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THE SOLUTION BY ITERATION OF NONLINEAR INTEGRAL EQUATIONS 


abstract 


In the mathematical formulation of many physical problems there are fre¬ 
quently encountered integral equations of the linear as well as non-linear 
type. The adequate solution of such equations then becomes a practical 
necessity. 

It turns out that for large classes of non-linear integral equations 
describing certain wave phenomena the well known method of successive approxi¬ 
mations provides a feasible means of constructing a solution, securing its 
uniqueness, and, in the process, furnishing practical estimates of relevant 
magnitudes and errors. Furthermore, the method also lends it3eir quite easily 
to treatment on high speed computing machinery. 

These principal conclusions are illustrated - and borne out - by means of 
an appropriately chosen example. 


1. INTRODUCTION 

It ia a well known fact that many physical problems lead to nonlinear integral 
equations. The following examples serve to illustrate this point. 

±. Nystrom has shown [lJ that the general first boundary value problem in the 
theory of ordinary differential equations may be reduced to an ordinary differential 
equation of the type 

y" - F 6c,y) -1/2 Jr x h 1/2 

with the boundary conditions y(-l/2) » y(l/2) ■ 0. This equation, in turn, is 
equivalent to the integral equation 

y 6c) - J K(x,a) F (s,y(s)) da, 

- 1/2 

X6c,a) representing the Green *3 function ">? the boundary value problem. t* 

I 

ii. boundary value problems of the form f 

" UBU xx + tt // " ex P u 

with u » f 6c,^ ) on the boundary C of the given region B i are reducible to integral 
equations, as follows [ 2]1 Let w(x, £) be the solution of the boundary value 
problem 

4v*0, w » f on C, 

Then the function y - u-w cleariy satisfies 

Ay-y exp w - [exp (y)-yj exp w 

with y - 0 on C. If now K(x,£ ; s,a)£ 0 denotes the Green's function of the problem 
4 y-y exp w » 0, y - 0 on C, 
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then it may be shown that, y(x, £ ) i 3 a solution of the nonlinear integral equation. 

y(x, £) * ff K - x > * > s f°) t r ) -yj ®:<P (w) dsd“ = 0, 

iii. In the discussion of certain wave phenomena L. Lichtenstein £3j was led 

;-,o the 3tudy of integral equations of the form 
„ b b 

L y (x) - X / / K., (s,s^.. .3 i >5(s^"y(s^) ds 1 ...ds i - 0, 

i»l a a 

with L denoting a linear operator, and the symmetric and continuous functions 
in the (n+l) dimensional unit cube. 

Physical situations leading to similar nonlinear integral equations have also 
be.'a described by Sc-hmeidler 0_i J - 

The solution of such integral equations may be accomplished in general by a 
variety ox methods, with the choice of the most appropriate one naturally depending 
upon the form and type of the integral equations of interest. Such methods may 
emply closed analytical expressions, power series expansions, the replacement of 
integrations by summations, successive approximations, variational methods, etc. J. 
The approach to be used here will be the iterative one, to be applied to 

integral equations of a generalised Lichtenstein types 
r b b 

y(x) - ^ (••• $ K i (x,s 1 ,...s i )F i (s 1 ,...S£) y(a 1 )...y(s i )'ds 1 ...ds i - 0. (l) 

i-1 *a a 

Integral equations containing an ’'inhomogeneous” term, e.g., of the form 

y* 00 - K q (x) ^ 2 * 2 (y*) * 0 

ean obviously be reduced to the form (l) by the substitution.- 

y - y - V 

For kernels whose squares are integrable, and functions F^ that have bounded 

derivatives the existence and uniqueness of a solution will be proved, and an 

upper bound for max |y| will be given. The proof is based on the method of successive 

iterations, a method which also easily provides an estimate of the error inherent in 

2 





the nth approximation. The results will then be Generalized to ay stems of equations 
of type (l), and, finally, an example will be given. 

9. ASSUMPTIONS .WO NOTATIONS. 

For the sake of simplicity, we shall take r = 2 in eq» (l), i«e. consider in 
the interval I: a k x,s,t h. b the equation, 

y(*)- J* 1 (x,s) F-,'s,y(s))ds -£fKg(x,s,t)F 2 (s,t,y(E»i - , ( t )) dsdt “ °* (2) 

The kernels and K 2 will be assumed to have integrable squares, so that J*K^(x,st)ds, 
(x, 3 »t)ds dt exist and are continuous functions of x £l. For the maximum of 
a continuous function f(s) defined in I ve shall employ the notation 
II ft - max |f 6c)/ . 

Xcl 

Then 

K fc * UfK^ da + Jj j K^ 2 ds dt 0 certaihly exists. 

The continuous functions F.(s,u), F-(«?,t,u,v), defined for s, t € I, u,v e , 
where J denotes a certain interval ^*c,c}, are assumed to satisfy conditions of uhe 
form 

I F,^ ) - F 1 (S^u 2 )|4 \ \\ “ «2 l 0) 

/F 2 (s,t,u 1 ,v 1 ) - F 2 (s,t,u 2 v 2 )j£ L 2 ( +K - v £ f 

where L-, Lq are positive numbers independent of s, t, u, v. The quantities 

Mj 2 “ II F 1 2 (a>0)/| , Mg 2 * i/F 2 2 (s, t, 0, 0) ji 

then also exist. 

Relations (3) certainly hold if, for example, F^ and Fg have continuous partiSL 
derivatives, 3F,/fcu, dFg/ou, bFg/bv for arguments restricted, respectively^ the 
internals I, J. 
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Further, we put 


A - b - a, G « 1 or 2 for A fz I or A >/ 

Other quantities to be employed later on are 

9 2 2 2 2 2 

L ~ - L 1 + h > M “ \ * m 2 -• 


A KM - E, A° KL - B. 


In the following extensive use will be made of Schwarz's inequality 

/J*f (x)g(x)dx| ± ( /f 2 dx) 1 / 2 ( /g 2 dx) l/2 . . (h) ' 

Thi3 relation is applicable also to multiple integrals and sums. If wo put 


Iff I m I^f 1 / 2 ); 

then eq. (U) may be expressed in the form 

i z (f 1/2 e 1/2 ) 4 Kf)i(g). (5) 


3. LEMMA. 

We are now in a position to sw.o the following Lemma: 
Let y 6J be ; 

(3). Then #y I (l - 2D)_£ E,._„ 

Proof. It follows from eq. (2) that 

'»*»■* Jfw 4 iff K 2 F 2 i' 


In this inequality the right side is,by (3), identically in x, 

4 10^) I (F^ ♦ I (Kg) 1 (F 2 ). 

From Schwarz's inequality for sums it may now be deduced that 
lyU ♦ I 8 (K 2 )J 1/2 [l 2 ',t\) * I S {F 2 )J 1/2 

* V 2 - 

Since, by assumption, y «- J, we may make use of relationships (3), getting 
I 2 (Fg) 4j^£/F 2 (s,t,o,0)( + L 2 (/y(s)/ + jy (t)j )J 2 dsdt 

#y« > 2 - 
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I 2 ^) k ^ 2u (H 1 + 2 L x ;/ y // )\ 

By ad?ition, then, 

i 2 ^) + I 2 ii 2 )& + 2L tty B ) 2 . 

Thus 

By Ilk K A°(K + 2Lfyf ) * E + 2Dgyg , 
which oroves the lemma. 

|j, EXISTENCE THSGREK. For functions F^, F 0 satisfying conditions (3) with 

0 *E/(l - 2D) * & (63 

eq. (2) has exactly one solution. 

To prove thi3 theorem we U3e successive iteration, starting with arbitrary con¬ 
tinuous function y Q fc J, and defining a sequence {^ of continuous functions by 

y n+ ivst) -jk,(x,s) F 1 (s,y n (s)) ds + JJ" K 2 (x,3,t)F 2 (s,t,y n (s), (7) 

y n (t) ds dt, n - 0,1,2... 

First we show by induction that y n € J. Let, then, y^_ n have this property. 
By definition (7), following the same reasoning that wa3 used to prove the lemma, 
there is obtained. 

tlyj £ E + ZD // y^f . 

From the inductive assumption and condition (6) it may now be inferred that 
Hy n fl*E + 2DOi: C, 

l.e. y n 

Next, the sequence is easily shown to be uniformly convergent, as follows: 

From 

i W^n I 4 - F l»n-l )i ds i 

* fiT' t 2 < F 2 !y n >‘ F 2 (y n-l )) d = dt / 

it follows, by means of Schwarz's inequality, that 


5 




|y n +l " : J~ r a-l J ' * 1 1 * 

Taking advantage of rci-t.i*iris (3) l'. - «ij<t si- 

- = k .*. a VVl //, 


so x^nax, 


or, by induction, 

90 

Since 1 - 2D >0 JJT (2D) is convergent, and it follows that the partial sums 


7 r ° f y o + ^ (y n “ y n-1 ) are absolutely 


and uniformly convergent. 


The continuous limit function y of the { Y n j is easily shown to be a solution 
of eq. (2). Putting 

J KjF (y) ds + Jf*2 F 2^ ds dt “ y *^ 
we find^namely, that 

I>Vl - /u it - ^W) * i 2 (F 2 (y n ) - f 2 (y)).] V2 . 

But the continuity of the in conjunction with the uniform convergence of y^ 

towards y suffice to make the difference, F^(y n ) - F^(y) arbitrarily small by 
taking n sufficiently large. Thus also lim y n “ y uniformly in z. It follows 
that y* • y, i.e. that y is a solution of eq. (2). 

Finally , to show that y is the only solution of eq. (2) under the conditions 
stated above let it be assumed that there exists another solution Y (::). Clearly 

Hw««*jr n+ 1 -r<l *«y n . 1 -yll , 
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and Hy^j-yH ^ c/2 for n rN^. The reasoning applied above to establish the 
relationship (o) no’-; leads to 

0 5V L -n*(n» K ty -ni, 

or 

l»7 c+i - ' r *±(™fHy 0 - y H> 

whence 

fl y n + l ” Y * 4: * /2 * or n >Ng« 

Consequently 

jf Y - y // < C . 

Due to the arbitrariness'of C and the meaning of |Y - y/fit is then necessary 
that Y(x) - y(x) * 0. 
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By the lemma # which if applicable here, and by condition (6), 
lj y//£ E/(l - 2D)j£ C+ 

For increasing values of K, L, M, therei‘orp,the upper bound 
E/(l - 2D) of Jyjalso increases. 

An estimate of ths error {/ y - y^jj of the n th approximation y R is 
easily obtained. Fran 

»7-y n n«wV *“ y -Vi" 

it follows that 

n+m 

n y - y„at 2 n*i.i'V - Wi u • 

i-n 


For fixed n. m-**>•. theny since y^-^y uniformly, 

lly-y„l± 1 «y ul - r ± S 

n-l 

■ - 2 " y i*i -yj' 

1 i-o 

l| 7 - yjtllyx - y J (2D) n /(d - 2D) for n£0. (9) 

5. GENERALIZATIONS. 


system 


Equation (2) may be considered to be a special case of a more general 


y(x ' ~ fj-Jh^’h . 


y P (8 j) ds i— dB j ■ °> 


where i - 1,2,..., p, p 2*1, q Z-l, To such systems the results of the 
previous sections may be transferred without difficulty; it is necessary 
only to replace the quantities o, E, D y l(, etc. by 





<i • 1 or q for A .61 or A > 1, 


I 



i! 1. i 2 ( K i<H 
J 


: K i M i - e, a“ l l 


*»• D, 


II y il - »/ ^ ^00 J * 


Finally, the important subsidiary condition 1 - 2D > 0 must be replaced by 1 ~qD >0* 
For p ■ q ■ 1, i.e. # for the equation 

y(x) *• fK(x,!i) p (s t y(s)) da » 0, (10) 


this condition I » qD > 0 becomes 0 AKL *: 1, or 

oa¥L 2 ^.i, (ll) 

Let us ascertain its meaning for the case of the linear integral equation 
of the second kind 

r-4~ 

y(x) - Xj K (x,s) y(s) d* - f (x). (12) . 

This type of equation is obtained from eq» (10* by putting 
F(a,u) ■ g(s) + X u 
f(x) «Jk(x,s) g (s) da. 

Since in this case L » II dF/du (I - (\ ( t inequality (11) indicates that 
k? ff K Z («*t) dadti lVk 2 t 1. 

Now for the characteristic value of smallest absolute amount of the 
symmetric kernel there exists the relation {6 J 

V * ^If . 

Consequently, condition ('ll.) becomes 

I X/. -c 1^1, 

which is the well-known sufficient condition for the existence of a unique 
solution of ea, (12) for arbitrary ±\x*» 

? 




6v APPLICATIONS. 

There are many practical cases where the kernels may be expressed in 
bilinear forms 

K-Jxjs) - (x '/'li (S) (W) 

¥ x » , ’ t) - f f 2} (x) f 2j (3.t) , 

the being linearly independent of the 0-... and the summation to be extended 
over a finite or infinite nunber of values i, j. In such cases the iterative 
procedure (7) clearly becomes 

* puM Vi.i *pt) <*> Vi.j '< 1W 

*W,i * /Vli <e) F i < s * % 4l (:5 >7V * V ) > d ‘’ 

V.j ‘ If Y?i * 3 (*»*» -f (•) ♦ J («). ^(t) ♦ at. 

The iterative procedure tor the functions 7 n (x) is thua reduced to a metaod 
of successive approximation for the constants p^, q^s 

Vl.t • *l< p nl' %3> < 15 > 

In the limit, then, 

y(x) « (*) P ± ♦ ZV 2 j (*) qj , (16) 

where the p it q^ satisfy the functional equations 

p i • f <p i> V 

for all i, J. 

q g - $ 2 (p ± , qj) (17} 
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Whenever ths functional equations (1Y) an.- directly solvable the 
solution (16) may be obtained immediately. Otherwise U3e must be made of 
the successive approximations (15)» 

Let us now apply the foregoing deductions to the following cases 
IC^XjS) « x(l-s), Is <o, 1^ 

K 2 (x,s,t) - exp(x-s-t) 

F^(Sj,u) ■ X(u+1) Js ^-1, 1/* 

F 2 (s,t,u,v) « X u v. 

The equation to be s olved, than, is ^ 

y(x) -X / /x(l-s) fy(s)+lj ds ♦ Jj”exp(x~ s-t)y(s)y(t)dsdt| - 0. (16) 

For the convergence of the Iterative procedure discussed above it is 
sufficient that 

0 £ E/(1-2D)^1. (6) 

Now 

A i 1/2 

K - £(1/3) ♦ (e 2 - l) 2 Ae 2 j * 1.31 
L « 2 1 / 2 J A | , M - / X I . 

Condition (6) is satisfied, then, for 

| X| * X q * [K(U2 |T5)J * 0.2 (19) 

further, the lemma indicates that for the solution y(x), 

/; y(y)fi= K ; xj /{I - 2K ,'Xj /F) . 

Since in this example i » j « 1, the iterative process for eq. (18) may be 
expleased in the fora 

y„l(x) - X <p n ,. Vx «cp (x>; 
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with 


P w i - J (1-s) [X(p n s+q n exp (s)) + \j da 

o *1 2 

q r»l ' x ‘-J <V K:ii <- a) * q n' ds J • 

0 

According to the existence and uniqueness theorem p - lira p_, 
q • lira q n exist. These limits satisfy a system of equations corresponding 
to equations (17). in this case 

p(l-X/6) “ qX(e-2) - 1/2 
X 2 £(1 -2/e)p « q] 2 - q. 

Eliminating p now there is obtained for q the equation 

x 2 (« p ♦ p) 2 = y q (20) 

a • 6X(«s-2) 2 > (6-X)e 

f * 3(«-2) 

JT- e 2 (6-X) 2 . 

Consequently, 

1/2 

q * ^(Y - 2aBX 2 ) ♦ fcfy - UapX 2 )J }>/ 2a 2 X 2 . (21) 

Of the two possible solutions indicated here the one corresponding to the 
positive sign becomes unbounded as X approaches zero. Since it is desired 
to keep the solutions for all / X/ h. X q of norm not exceeding unity the positive 
sign in eq. (21) must be discarded. 
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Suppose we wish to determine the solution of eq. (18) for ). <• 0.2. 
For the initial approximation we take 


y o (x) - 0.2 (x + exp (x))* 
The raaiilting approximations 

Vl • °- S r p n /d 


• 1/2 


Vi ■ 0 - 0il * %J 

are shown in the table below. 

“Vole I. Successive Approximation 


n 


q n 

0 

1.0000 

1.0000 

1 

0.6770 

0.0639 

2 

0.531? 

C.CC2«* 

3 

o ,5 i«i 

0.0008 

h 

0.5171* 

0.0008 


The exact solution, as determined from eq. (21), is 

y(x) - 0.2 (0.51735 X * 0.00075 exp (x)). 

While there is always convergence to the solution for y Q € J, IXI i- \ Q , 
the choice of an initial approximation y Q outside of J, or a value of J XJ 
say lead to divergence. This is the case, for example, for p Q ■ q Q ■ 25, 
X * 0.2, and p - q - 1., X - 1„ 

MASK LOTKIN 
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It is worthwhile mentioning that iterative procedures such as those 
exemplified by eqs. ( 15 ) are eminently suited for treatment on high speed 
computing machinery. The required quadratures may be performed uy anj 
suitable standard methods, or by means of the formulas outlined in f?JL 
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